driver_module Module


Uses


Variables

Type Visibility Attributes Name Initial
integer, public, parameter :: length = 5000000

common/agh/xnpar,vel,dland,dcoll,dalf,perpn,tetai,jrad,iww,izz
integer, public :: irs
integer, public :: iabsorp
real(kind=wp), public :: rzz

common /abc/ rzz,tetzz,xmzz

real(kind=wp), public :: tetzz

common /abc/ rzz,tetzz,xmzz

real(kind=wp), public :: xmzz

common /abc/ rzz,tetzz,xmzz

integer, public :: iznzz

common /abc/ iznzz,iwzz,irszz
integer, public :: iwzz

common /abc/ iznzz,iwzz,irszz
integer, public :: irszz

common /abc/ iznzz,iwzz,irszz
real(kind=wp), public :: hrad
integer, public :: ind
integer, public :: im4
real(kind=wp), public :: pow

common /acg/ pow
real(kind=wp), public :: pintld4
real(kind=wp), public :: pintcl4
real(kind=wp), public :: pintal4

Functions

public function dql1(ifound, jfoundr, pabs) result(powccc)

pintld=-pintld4*dfdv

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: ifound
integer, intent(in) :: jfoundr
real(kind=wp), intent(in) :: pabs

Return Value real(kind=wp)


Subroutines

public subroutine memorize_trajectory_point4(ro, theta)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: ro
real(kind=wp), intent(in) :: theta

public subroutine memorize_trajectory_point(vz, j, ro, powccc, driver)

memorize trajectory point

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: vz
integer, intent(in) :: j
real(kind=wp), intent(in) :: ro
real(kind=wp), intent(in) :: powccc
integer, intent(in) :: driver

public subroutine driver2(ystart, x1, x2, xsav, hmin, h1, pabs)

solve eqs. starting from xbeg ystart(1) = tet ystart(2) = xm x1 = xbeg rini x2 = xend

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(inout) :: ystart(2)
real(kind=wp), intent(in) :: x1
real(kind=wp), intent(inout) :: x2
real(kind=wp), intent(inout) :: xsav
real(kind=wp), intent(in) :: hmin
real(kind=wp), intent(in) :: h1
real(kind=wp), intent(in) :: pabs

public subroutine extd2(x, y, dydx)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x
real(kind=wp), intent(in) :: y(:)
real(kind=wp), intent(inout) :: dydx(:)

public subroutine difeq(y, dydx, nv, x, htry, eps, yscal, hdid, hnext, derivs)

! зачем save ?????

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(inout) :: y(nv)
real(kind=wp), intent(in) :: dydx(nv)
integer, intent(in) :: nv
real(kind=wp), intent(inout) :: x
real(kind=wp), intent(in) :: htry
real(kind=wp), intent(in) :: eps
real(kind=wp), intent(inout) :: yscal(nv)
real(kind=wp), intent(inout) :: hdid
real(kind=wp), intent(inout) :: hnext
procedure(Iderivs_func) :: derivs

public subroutine mmid(y, dydx, nvar, xs, htot, nstep, yout, derivs)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: y(nvar)
real(kind=wp), intent(in) :: dydx(nvar)
integer, intent(in) :: nvar
real(kind=wp), intent(in) :: xs
real(kind=wp), intent(in) :: htot
integer, intent(in) :: nstep
real(kind=wp), intent(inout) :: yout(nvar)
procedure(Iderivs_func) :: derivs

public subroutine rzextr(iest, xest, yest, yz, dy, nv)

rational extrapolation зачем save ????

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: iest
real(kind=wp), intent(in) :: xest
real(kind=wp), intent(in) :: yest(nv)
real(kind=wp), intent(inout) :: yz(nv)
real(kind=wp), intent(inout) :: dy(nv)
integer, intent(in) :: nv

public subroutine pzextr(iest, xest, yest, yz, dy, nv)

polynomial extrapolation

Arguments

Type IntentOptional Attributes Name
integer :: iest
double precision :: xest
double precision :: yest(nv)
double precision :: yz(nv)
double precision :: dy(nv)
integer :: nv

public subroutine driver4(ystart, x1, x2, rexi, hmin, derivs)

var if(dabs(dydx(3)).ne.zero) h=dabs(hmin1/dydx(3))/hdrob1 if(dabs(y(3)-rexi).gt.rrange.or.nstp.eq.maxstep4) then ! exit !sav2008

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(inout) :: ystart(:)
real(kind=wp), intent(inout) :: x1
real(kind=wp), intent(inout) :: x2
real(kind=wp), intent(in) :: rexi
real(kind=wp), intent(in) :: hmin
procedure(Iderivs_func) :: derivs